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Perturbative and non-perturbative QCD 
C.J. Maxwell 

Institute for Particle Physics Phenomenology (IPPP), Durham University U.K. 



In these lectures we give a concise introduction to the ideas of renormalon calculus in QED and QCD. We 
focus in particular on the example of the Adler D function of vacuum polarization, and on relations between 
perturbative renormalon ambiguities and corresponding non-perturbative Operator Product Expansion (OPE) 
ambiguities. Recent work on infrared freezing of Euclidean observables is also discussed, 
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1. Introduction 

In these lectures I will discuss the large-order be- 
haviour of perturbation theory. The field theory per- 
turbation series in the renormalised coupling is not 
convergent, and in nth order the coefficients exhibit 
n! growth. The sources of this divergent behaviour 
are instantons, which are connected with the combi- 
natoric growth in the number of Feynman diagrams, 
and renormalons which are associated with individual 
diagrams containing chains of fermion bubbles. We 
will focus on vacuum polarization in QED and QCD, 
and will analyse the Borel plane singularities which 
arise due to these chain diagrams, so-called ultraviolet 
(UV) and infrared (IR) renormalons. In QED the UV 
renormalons render perturbative results ambiguous, 
whereas the ambiguity in perturbative QCD arises 
from IR renormalons. The IR ambiguities in pertur- 
bative QCD are associated with non-logarithmic UV 
divergences in the non-perturbative Operator Product 
Expansion (OPE), and can cancel between the per- 
turbative and non-perturbative sectors. We shall par- 
ticularly focus on the deep connections between the 
perturbative and non-perturbative sectors, reporting 
in the final section on some recent work on infrared 
frezing [1,] which suggests that IR and UV renormalons 
conspire between them to provide finiteness and conti- 
nuity avoiding a Landau Pole in the coupling. For an 
excellent review on renormalon calculus see Ref.Q, 
and for a review on connections between perturba- 
tion theory and non-perturbative power corrections 
see Ref . 0] . Recommended texts on quantum field the- 
ory are Refs.fJ, 's', We begin with a brief introduc- 
tion to QCD, and introduce the vacuum polarization 
on which we will largely focus. 



2. Introduction to QCD 

Quantum Chromodynamics (QCD) is a non-abelian 
gauge theory of interacting quarks and gluons. The 
gauge group is SU (iVc), and there are — 1 gluons. 
Experimental indications are that Nc = 3. The La- 
grangian density is 



QCD 



Here a ~ 1, 2, 3, . . . , 8, and Tj, — Aa/2 are the genera- 
tors of SU(3), the 3x3 Gell-Mann A-matrices, satis- 
fying 



[Ta , Jfe] = ifabcTc 



(2) 



The quark fields carry colour, R, G, B, and transform 
as a triplet in the fundamental representation 



(3) 



C-QCD is invariant under local SU(3) gauge transfor- 
mations 



i^ix) — > U{x)'ip 



Ax) 



(4) 



The field strength tensor G^^ contains the abelian 
(QED) result and an extra term proportional to the 
structure constants f"''"^ which are responsible for 
three and four-point self-interactions of gluons, not 
present for photons in QED. 



abc/-ib 



(5) 



For QCD (but not QED) one also needs to include 
unphysical ghost particles. These are scalar Grass- 
mann (anti-commuting) fields needed to cancel un- 
physical polarization states for the gluons. The re- 
quired Fadeev-Popov extra term in Cqcd is 



C 



ghost 



(6) 



In both QED and QCD one needs also to include a 
gauge fixing term if inverse propagators are to be de- 
fined. 

1 



gauge— fixing 



2? 



(7) 



There is only one other gauge-invariant structure that 
we could add involving the dual field strength tensor 



Co 



647r2 



(8) 



This is a total derivative and so produces no ef- 
fects at the perturbative level. However, ii 9 
non-perturbative effects would induce a CP-violating 
electric dipole moment for the neutron, experimental 
constraints on this provide a bound 16*1 < 3.10^^°. 
■(1) 
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2.1. QED vacuum polarization 

A crucial ingredient in our later discussions will be 
the one-loop vacuum polarization diagram shown in 
Fig.l. Using the Feynman rules the diagram is given 
by 



iW{k'^) = (-1) 



J/3/3' 



[{k +py — nn?] 



(9) 



Here e is the QED coupling, a — e^/47r. The integral 
is divergent and requires "regularization" . The most 
widely applied method is "dimensional regularization" 
, in which the integral is performed in D = A — 2e 
spacetime dimensions and then the limit e — > is 
taken. We need to evaluate 



{2tt)D p^k+p)^ ' 



(10) 



The continuation to ZJ-dimensions endows the dimen- 
sionless coupling e with a mass dimension, [e] = 



2 — -J = £, so one needs to replace e 
finally obtains 



iW^ik^ 



Stt' 



■j''''k^] 



In- 



- (ln47r-7£;) - - 



One 
(11) 



Here the divergence as e — > is contributed by 
r(e) = i _ + 0{e) , where je ~ 0.5722 is the 
Euler constant. Counterterms are introduced to re- 
move the 1/e divergences (see next section) , the finite 
contribution they also cancel is arbitrary and deter- 
mines the subtraction procedure. Modified minimal 
subtraction (AIS) absorbs the ln47r — 7^ term, mini- 
mal subtraction (MS) does not. 



to 1/e, the divergent parts present in loop calculations 
can be cancelled. One finds that the renormalised 
coupling runs logarithmically with the renormaliza- 
tion scale /i, and satisfies the beta-function equation 



da 
din 11 



/3(a) = -6a^(l + ca + C2a^ + csa^ + . . .)(13) 



Here a(^^) = as(^^)/7r ~ 5s(/i^)/47r^. The terms up 
to and including C3 have been computed in the MS 
renormalization scheme. The first two coefficients are 
scheme- independent 



b = 



{UN, - 2Nf) 



(153-mV) 

126 ^ ' 



The first TVc-dependent contribution to b arises from 
gluon and ghost vacuum polarization contributions, 
and the second iVy-dependent contribution is (up to 
a group theory factor of T{R) ~ 1/2) just the QED 
vacuum polarization contribution considered earlier. 
Nf is the number of active quark flavours (fermion 
species in QED). For SU(3) QCD with Nj < 33/2, b > 
0, and a{fi^) ^ as /^^ — > cx) (Asymptotic Freedom). 
At the one-loop level the solution of the equation is 



61n(/i2/A2)' 



(15) 



At two-loops the solution may be written in terms 
of the Lambert W^-function f7l| defined implicitly by 

Wiz)exp{W{z)) = z 



1 



c[l + w{z{^i)] 



- -US'"' 



(16) 



At higher loops a{fi^) will depend on the choices of 
the non- universal beta- function coefficients 02,03,.... 
A special choice is an Hooft scheme [i] where C2 = 
C3 . . . = c„ = in which case a^fi^) may be written in 
terms of W{z) as above. 



2.2. Renormalization and running 
coupling 

One needs to distinguish between the bare La- 
grangian Cb and the renormalised Lagrangian Cr. 
Bare and renormalized fields and couplings are related 
by infinite renormalization factors Zi. In massless 
QCD for instance one has 

TpB ^ zl^^-tJjR , G^'" = Zg^^G'j^'" ,gsB ^ ZggsR.{l2) 

This infinite reparametrization can be implemented 
by introducing counterterms into the Lagrangian. 
These will contain counterterm coefficients, by choos- 
ing suitable values for these coefficients proportional 



2.3. Vacuum polarization and the e+e~ 
hadronic total cross section 

Consider e+(pi) -I- e~(p2) X, where A" is a 
hadronic system having total momentum q = pi + P2- 
The squared e~'"e~ cm energy is s = q^. The relevant 
amplitude is 

27re2 

A{e+e- ^ A) = v{p,)ru{p2){X\J^\Q). (17) 

s 

Here is the electromagnetic current for quarks 

J^ = J2Qf^fr^f, (18) 
/ 
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p + k 

Figure 1: One-loop diagram for QED vacuum polarization 



summed over active quark flavours. The e'^e total 
hadronic cross section is then given by 

a ^ ^a{e+e- X) 

X 

* X 

X (X|J''(0)|0)(X|J'^(0)|0)*. (19) 
Here the leptonic tensor l^,^ is given by 



(20) 



Using the optical theorem the term may be related 
to the imaginary part of the photon propagator. We 
may write 

Ex {2nts(^\p, +P2- q){0\MO)\X){Q\JAQ)\Xr 
dW«-^(0|[J^(a;), J,(0)]|0) 

= i^-n-(.). 

Here H^''(Q^) is the vacuum polarization function 



(21) 



Uf-^iQ^) ^IGn^i J d^xe"J-^{0\T{J^{x)J„{Q))\0){22) 

Here = — > 0. Conservation of J^, S^J^ = 
then dictates the tensor structure 



(23) 



The quark-loop vacuum polarization diagram we con- 
sidered earlier contains hi{—s/fj?) which has an Im 
part of ITT. Combining with QCD colour factors, and 
taking into account the fractional charges of quarks 
one finds 



47r, 



(24) 



/ 



One then finds for the total hadronic cross-section 



a{e^e — > hadrons, s) 



3s 



-3J2Q}- (25) 



A convenient observable to measure at e~^e colliders 
such as LEP is the Re+e- ratio, defined by 



(T(e+e —t hadrons, s) 
(T(e+e^ s) 



(26) 



The /i^/i cross-section can be directly measured in 
the experiment, 



cr e e 



^l+p ,s) 



3s 



(27) 



The total hadronic cross-section differs from the /i^/i^ 
point cross-section by the factor 3 (= Nc colours for 
each quark/anti-quark), and which takes into 

account the fractional quark charges. Taking the ratio 
one finds 



f 



(28) 



Here TZ{s) = TZpt{s) + TZnp{s) corresponds to QCD 
corrections to the parton model result. The perturba- 
tive (PT) contribution is of the form 



UpTis) = a{s) + ria^(s) + r2a^{s) + 



(29) 



Here a(s) denotes the coupling in the MS scheme, 
the ri and r2 corrections have been computed in the 
MS scheme. The non-perturbative (NP) component 
arises from the Operator Product Expansion (OPE). 

Only n((5^) — H(0) is observable and so it is conve- 
nient to take a logarithmic derivative with respect to 
and define the Adler D function. 



(30) 



This has the same parton model expression as 
Re+e-{s) with the perturbative corrections TZpt{s) 
replaced by VpTis), 

VpTiQ^) = a{Q^) + dia\Q^) + d2a\Q^) + ... .(31) 

7?.(s) is related to I?(— s) , by analytical continuation 
from Euclidean to Minkowskian momenta. One can 
write the dispersion relation 



n{.s) 



1 , V{t) 



(32) 



IPM-LHP06-sch 



4 



IPM-LPH06, Tehran, Iran, May 2006 



One can also write this as a circular contour integral 
in the complex-s plane 



deV{se''^). 



(33) 



3. Large-order behaviour of PT- 
Instantons and Renormalons 

In 1952 Dyson presented an argument [§, that 
QED perturbation series are divergent series with co- 
efficients growing like n\ in nth order. Consider 



fie') = E/"^ 



(34) 



Assuming that /(e^) converges for > 0, implies 
that /(e^) is analytic at = 0, and must converge 
for some negative values of e^. For < 0, however, 
like charges attract and unlike charges repel ! The 
vacuum state is then unstable and it becomes ener- 
getically favourable to create more and more e^e^ 
pairs. Consider a system of N interacting electrons, 
the energy will be 



E ^ NT + ^N^Ve^ 



(35) 



Here T is the mean kinetic energy, V the mean 
coulomb potential and ■^{N'^ — N) ~ ^TV^ counts the 
interacting e~'"e~ pairs. For > the energy (wrt 
N) is bounded from below and there exists a stable 
minimum, the vacuum state, at TV = 0. For < 0, 
however, there is an initial minimum at = be- 
yond which E rises ^ N until a maximum is reached 
atN = N„ru 



Nr. 



T 



1 



V\e\ 



(36) 



Beyond this point E decreases as — A^. There is no 
stable minimum. One infers that the divergent nature 
of the perturbation scries emerges when more than A^c 



terms are considered, 
decrease and /„e^" ~ 



For n < Ncrit the /ne terms 



fn+ie 



2n+2 



at n '--^ Nc- 



fn+l 1 



fn ^ ni 



(37) 



3.1. Asymptotic Series/Borel summation 

Consider a function f{g) expressed in terms of a 
power series expansion 



fig) 



n=0 



(38) 



Consider a domain T> of the complex g-plane such that 
arg g < tt/2. The series is said to be asymptotic inside 
T> if the series diverges for all g 7^ and 



N 



/(5)-E/"f" 



n=0 



< fN+l\9\ 



N+1 



(39) 



The crucial property of asymptotic series is that the 
error made in truncating the series is less than the 
first neglected term. If /„ ~ n! then successive terms 
will decrease until a minimum is reached at Nopt 
terms, thereafter the size of the terms will increase 
without limit. By truncating the series at rt = Nopt 
the best possible approximation is found. 

One can define a function which is asymptotic to 
the series by using the Borel method. If /„ ^ nl one 
defines the Borel transform of the series 



mi^)-T.^9" 



(40) 



This series will now have a finite radius of convergence. 
One can then write 



fig) = / dze-^'/a E 



— z" 
n! 



(41) 



This follows since if one performs the integral term by 
term and uses the result 



dze-^/^z" = nig". 



(42) 



one formally reproduces the divergent power series for 
fig). If this series has a finite radius of convergence 
then one can show that the B[f]{z) series has infinite 
radius of convergence , and f{g) is equal to the Borel 
sum 



/•OO 

/(g)- / dze-^/^B[f]{z). 

J O 



(43) 



If f{g) has zero radius of convergence with /„ ^ nl 
then B[f]{z) will have finite radius of convergence, 
and can be analytically continued outside of that ra- 
dius to the whole of the integration range [0, 00] 



fig) 



dze-^/3B[f]{z). 



(44) 



Here the « symbol means "is asymptotic to" . Notice 
that there is in general not a unique function to which 
the series is asymptotic since we can always add a 
term e~'-^/^, which has an identically zero Taylor 
expansion in powers of g, and the series will also be 
asymptotic to that function. If we have information 
on the analytic structure of f{g) in the complex 
g-plane it is sometimes possible to enforce that C = 
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(Watson's Theorem). 

Two relevant examples which will be important 
later when we introduce renormalons are alternating 
and fixed-sign factorial growth. First consider the al- 
ternating factorial series, 



(45) 



Tl = 



This has the Borel transform 

oo 

S[/-](^) = E(-1)' 



(46) 



n=0 



This series converges for \z\ < 1, and may be analyti- 
cally continued to B[f^]{z) = 1/(1 + z) on the whole 
range [0, oo] one then finds, 

f°° 1 
/_(z)«/ dze-^/s-— ^ = -Ei(-l/5). (47) 
Jo (1 + 2) 

Here Ei(x) is the Exponential integral function defined 
(for a; < 0) as 



Ei(a;) = - / dt 

J-x ^ 

Now consider the factorial series 

oo 
n=0 

This has the Borel transform 

oo 

B[/+](^) = ^z". 

n=0 



(48) 



(49) 



(50) 



This series converges for \z\ < 1, and may be analyti- 
cally continued to B[f^]{z) = 1/(1 — z), 

/+(<?)-/ dze-^/s- r-Ei(l/g). (51) 

Jo (1 - 2) 

For a; > Ei(2;) is defined by 

poo t 

Ei(a;) = -PV / dt . (52) 

J-x ^ 

There is a pole at z — 1 in the Borel plane which 
renders the result ambiguous, with a ±i7r contribution 
depending on whether the integration contour passes 
above or below the pole. Use of the principal value 
(PV) prescription is equivalent to averaging over these 
choices. 



3.2. Proliferation of Feynman diagrams 
and Instantons 

Consider the following integral expressed in terms 
of a power series expansion 



I{9) = 



1 



+ 00 



dx e-(-V2+9.V4) 



(53) 



n=0 



This integral is the generating function for the number 
of Feynman diagrams contributing to the vacuum-to- 
vacuum transition amplitude of 0"* field theory. One 
has 



In = 



+ 00 



j_{-ir 

^/27^ 4"n! 

(-l)"r(2n+l/2) 
r(n+l) 

(_4)" 



2tt 



-{n-iy. 



(54) 



and so the combinatoric growth in the number of 
Feynman Diagrams would be expected, other things 
being equal, to contribute to factorial growth of the 
coefficients. 

We now turn to a brief discussion of instantons [l3] ■ 
Consider a generic Green function Q (g) for a simplistic 
field theory of a single field at a single spacetime point. 



1 f+°° 

Gig) - - / ' 
g J-oo 



d(j) e 



-S{4>)/g 



(55) 



This expression can be written in terms of the Borel 
transform of Q{g). 



Gig) 



1 

9 Jo 



-"/f B[g]iz)dz. 



(56) 



By inspection, the Borel transform of Gig) is found to 
be 



B[g]{z) 



+ 00 



• siz~sm, 



(57) 



and this can be rewritten by change of integration 
variable as. 



/ + 00 
d5(0) 
oo 



dSi(b) 



6{z - S{cb)) 



E 





-1 


34) 





(58) 
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Here (t>i label all solutions (f) obeying S{(t>) = z. We 
can see that singularities in the Borel plane occur at 
values of S'[0] for which satisfies 



0. 



(59) 



Hence (pi represent extrema of the action and they 
are therefore solutions to the classical equations of 
motion. Consequently, singularities exist in the Borel 
plane at values of the action corresponding to these 
classical solutions. We expand the action in the region 

of 

S{4>) - S{<|)i) + ]^S"{<|y,){<j,-(t>^f. (60) 
We then rearrange this whilst defining Zi — S{(j)i). 

1/2 



55(0) 



, (5(0) -5(00) 

' S"{<|>^) 
1 1 



Substituting this result back into the Borel integral of 
Eq.(56) one finds 



n=0 



Gig) = 

e-^/3(l-z/z,)"^/^rf2- (62) 

The series coefficients are then determined to be 
l\ r(n+ 1/2) 



^^^J r(i/2) ' 

' \ n 



1 

Zi 



n 



1/2- 



(63) 



Thus one finds nl growth in perturbative coefficients. 
The Borel plane singularities are at positions in the 
z-plane corresponding to the values of the action rep- 
resenting instanton solutions. By considering specific 
examples such as the anharmonic oscillator [lO| one 
can make a concrete connection between this n! in- 
stanton contribution and that due to the proliferation 
of Feynman diagrams discussed above. 

4. Large- A^/ approximation for vacuum 
polarization 

Consider the Adler _D-function we discussed ear- 
lier (Eq.(30)) with perturbative expansion. The nth 
perturbative coefficient coefficient d„ in Eq.(31) may 



be expanded in powers of Nf the number of quark 
flavours 

d„ = 4"liVj^ + Ar«-i + . . . + . (64) 

The leading large- iV/ coefficient d„ may be evalu- 
ated to all-orders since it derives from a restricted set 
of diagrams obtained by inserting a chain of fermion 
bubbles inside the quark loop (see Fig. 2). A crucial 
ingredient is the chain of n-bubbles, B'^^^{k'^). This 

may be defined as a product of bubbles Hp^(/c^) and 



propagators Pfj.u{k 



-iP„ 



.{g^Liy - fc^(l -C)k,j.kv) 



(65) 



- ' ^2 

Stringing the bubbles and propagators together gives 



k=l 



-iP-^^+^^im) 



(61) Here ji = ax. To evaluate this we will need the results 



and 



n 



H° 



k=l 

We finally obtain 



(67) 



(68) 



(-i)"(no)' 



k^ 



kt^k" 
fc2 



B^;:^{e) = (-f)"(^o)"HP'^'^(A:^e = o)]. (69) 

This is proportional to the Landau gauge (^ = 0) 
propagator. Explicitly 



j^fiiy _ {k^g^^ kfj_k^) 



[kr 



rh4 



C 



(70) 



The constant C depends on the subtraction procedure 
used to renormalise the bubble. With MS subtraction 
(7 = — |. We shall choose to work in the "V-scheme" 
which corresponds to AIS with the scale choice /i^ = 
e^^^^Q^, in which case C — 0. Applying the Feynman 
rules to the three diagrams then gives d!r''a"+^ 

d*k d^p 
(27r)4 (27r)4 
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Figure 2: Bubble-chain diagrams contributing to d 



["] 



1 



(71) 



The loop integrals can be evaluated using the Gegen- 
bau er p olynomial x-space technique, with the result 

mm 



-2n - 



n + 6 

2"+2 



16 



(1-2 



+ l>m>0 
2m-n-2 



)C2m+l n\ 



(72) 



4.1. Leading-^ approximation and QCD 
renormalons 

The large- A^/ result of the last section can describe 
QED vacuum polarization, but for QCD the correc- 
tions to the gluon propagator involve gluon and ghost 
loops, and are gauge (^)-dependent. The result for 
no(fc2) is proportional to ~Nf/3 which is the first 
QED beta- function coefficient, b. In QCD one expects 
large-order behaviour of the form dn ^ KrC {h/2)^n\ 
involving the QCD beta-function coefficient h = 
(33 — 2Nf)/6, it is then natural to relace Nj by 
(33/2 — 36) to obtain an expansion in powers of b 



(73) 



The leading-6 term dn 



(L) 



-3)"dir'6" can 
then be used to approximate dn to all-orders, and an 
all-orders resummation of these results performed to 
obtain 'D^p^{Q^). In order for the result to be renor- 
malisation scheme independent we need to use the 



one-loop form of the couphng a{Q^) = 2 / b\n{Q'^ / A'^ , 
and in what follows we will work in the V-scheme dis- 
cussed earlier so that C = in Eq.(70). If we use 
the Borel method to define the all-orders perturbative 
result we obtain 

/•OO 

V^^^iQ^)^ dze''/''('^"^B[V^^^]{z) . (74) 
Jo 

The Borel transform is given by [l3| 

n\T,(^h( \ - ^ Mn)-Ai{n)zn , Ai(n)z„ 
n[Lypj.\(Z) — > -f- 

(l + t) + t) 

Bo{n) + Bi{n)zn Bi(n)z 



E 



n=l 



The residues are given by 

(-l)"+i(3n2 + 6n + 2) 



1 - ^ 



7.5) 



Ao{n) 
Aiin) 



3 n2(n + l)2(n + 2)2 

8 b(-l)"+H» + |) 
3n2(n-M)2(n-h2)2 



So(l)-0, Bo(2)-l, Bo(n) = -Ao(-n) n>3 
Si(l)=0, Bi(2)--^, Si(n) - -Ai(-n) n 0S) 

In Fig. 3 we show the Borel plane singularities. For 
the Adler function in leading-^ approximation there 

are single and double poles in B\D^p^]{z) at positions 
z = Zn and z = —Zn with Zn = 2n/h, n = 1,2,3. 
The singularities on the positive real semi-axis are the 
infrared renormalons, IRn and those on the negative 
real semi-axis are ultraviolet renormalons, UVn- We 
shall see that they correspond to integration over the 
bubble-chain momentum fc2 in the regions < 



IPM-LHP06-sch 



8 



IPM-LPH06, Tehran, Iran, May 2006 



i 


Im z 


Z 


UV Renormalons 

A 




IR Renormalons 

A 


r ^ 

-8 -6 -4 -2 
' ' ' h h h h 


r 

4 6 8 
h b b 


... Re 2; 


• t • •! • • r • r • t • •! • • r • r • 


-| - - 1 - - |- • ^ - -| • • 1 • - |- - ^ - -| 

V J 

Instantons 



Figure 3: Borel plane singularities for the Adler D-function. Nf — 5 active quark flavours are assumed. 



and > Q^, respectively. The IRi singularity is 
absent for T> for reasons which we shall shortly discuss. 
Hence the singularity nearest the origin is the UVi 
renormalon, which generates the leading asymptotic 
behaviour 



27 

The IRn singularities lie on the integration contour, 
and hence there is an amiguous imaginary part de- 
pending on whether the contour is routed above or 
below the pole. This is structurally the same as terms 
in the non-perturbative Operator Product Expansion 
(OPE). We see that there are also singularities corre- 
sponding to instanton contributions. These are at po- 
sitions Zn = 4n corresponding to the nil pair contri- 
bution. The leading // instanton singularity at z = 4 
lies well to the right of the leading IR2 singularity and 
so does not dominate the asymptotic behaviour. 

4.2. OPE and IR renormalon ambiguities 

The regular OPE is a sum over the contributions 
of condensates with different mass dimensions. In the 
case of the Adler function the dimension four gluon 
condensate is the leading contribution 



Go(a(Q')) 



1 



(0|GG|0)CGG(a(Q')) , (77) 



where CaaiaiQ'^)) is the Wilson coefficient. The OPE 
is of the form 



(78) 



The n^^ term in this expansion will have the structure 
Cn{a{Q^)) - C„[a(g2)]'"(l + 0{a)) . (79) 



The exponent 5„ corresponding to the anomalous di- 
mension of the condensate operator concerned. Non- 
logarithmic UV divergences lead to an ambiguous 
imaginary part in the coefficient so that Cn = 
± idn\ If one considers an /i?„ renormalon 
singularity in the Borel plane to be of the form 
i^„/(l — zj Zn)''^ then one finds an ambiguous imagi- 
nary part arising of the form 



r(7n) 



g-.,./a(Q^)^l- 



1l + O(a)](80) 



Here the ± ambiguity comes from routing the contour 
above or below the real z-axis in the Borel plane. This 
is structurally the same as the ambiguous OPE term 
above, and if Cn^ = i^„7rz^"/r(7„) and (5„ = 1 — 7„, 
then the PT Borel and NP OPE ambiguities can po- 
tentially cancel against each other Taking a PV 
of the Borel integral corresponds to averaging over 
the ± possibilities. Notice that there is no conden- 
sate of dimension two in the OPE, the leading gluon 
condensate being of dimension four. This explains the 
absence of the IRi singularity at z — 2/6 evident in 
Fig.3. 



5. Freezing of Euclidean 
observables/skeleton expansion 

We wish to investigate the Q^-dependence of Eu- 
clidean observables such as the Adler function, and in 

particular whether V<^^\Q'^) = V^p^iQ^) + V^^Uq^) 
can remain finite as 0. We shall find that these 

components in the one-chain QCD skeleton expansion 
both vanish as ^ [1]. Let us introduce two 
other observables. The polarised Bjorken (pBj) and 
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Figure 4: Large- A'^/ contributions to tiie DIS sum rules at nth order in perturbation theory 



GLS sum rules are defined as [H, [l^ 



pBj 



97 ''"{x,Q'^)dx 



chain in a tree-level diagram, rather than a one-loop 
one. There are a finite number of simple poles. Using 
the integrals 



9A 



9v 



GLS 



1 - -CfJC{Q^) , 



(81) 



oo p~z/a 



(1 + z/'^n) 



= -z„e^"/'^Ei(-z„/a) 



oo p-z/a 



1 + — e-^"/'^Ei(z„/a) 



,(86) 



(82) 



/o (l + ^/^n)2 

we can obtain the following resummed expressions p!3 | 

oo 



K.{Q^) being the QCD corrections to the parton model 
result. We have neglected contributions due to "light- 
by-light" diagrams which when omitted render the 
perturbative corrections to Kqls and KpBj identi- 
cal. Finally, the unpolarised Bjorken sum rule (uBj) 
is defined as dOl 



1 - -CpUiQ^) . (83) 



Leading-^ resuhs for J^p^Q"^) and UJ^iQ'^) can be IC^p^Q^) = — 



computed from the diagrams in Fig. 4. The expres- 
sions for B[K}pl]{z) and B[u''p^]{z) are [2ll,[22 



4/9 



1/18 



8/9 



B[U^P^]{z) 



1 + - 
5/18 

Fi) 

1/6 

1 + i 

1/2 
1 - ^ 



1 + - 

Z2 



4/3 



1 - ^ 

21 



1 - ^ 

21 



(84) 



(85) 



These expressions are significantly simpler than for 
the Adler function since one is inserting the bubble 



a(Q2) 



^(Ao(r7,) - ZnAi{n)) - 2„^i(n) 



a(Q2^ 

-f (^o(?l) - Zn^l(n))] 

OO 



^{Bn{n) + ZnBi[n)) - z„Bi(n) 

(87) 



(Bo(n) + z„Bi(n))] 



1 

96 



Ei(-zi/a(Q^)) 



2e^^/'^(^')Ei(-Z2/a(Q^)) 
16e--i/'^(Q')Ei(zi/a(Q2)) 

10e--2/'^(Q')Ei(z2/a(Q2)) 



(88) 



1 

36 



8e-^i/'^(Q )Ei(zi/a(Q2)) 

6e-"^/"WEi(z2/a(Q')) 
2e^^/'^(«')Ei(-Z2/a(g')) . 



(89) 



These expressions have the property that they are 
finite and continuous at ~ where the one-loop 
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coupling a{Q^) = 2/b\n{Q^ /A^) has a Landau pole. 
For < the standard Borel representation breaks 
down and as we shall argue needs to be replaced by 
the modified Borel representation 



dze-''''^^"^B[vf^]{z) . (90) 



We shall show that this pair of Borel representations 
are equivalent to the one-chain skeleton expansion 
term in the two regions. The < representa- 
tion has an ambiguous imaginary part due to the UVn 
singularities on the integration contour. In Fig. 5 we 
show the behaviour of the three Euclidean observ- 
ables. The perturbative corrections are seen to change 
sign in the vicinity of = A^, and there is then a 
smooth freezing behaviour and they vanish as ~> 0. 
The finiteness is delicate. For small x 



Ei(2;) ^ln\x\+-fE + 0{x) 



(91) 



and so as A^ there is potentially a In a diver- 

ence. For V the coefficient of this divergent term is 
il 



J2zl[A,{n)+B,{n)] 



(92) 



n+l 



For /C^((5^) and Z-/p^^ (Q^) the equivalent coefficients 
are (-8 2 = 16 - 10 = 0) and (8 - 6 - 2) = 0, 



respectively. There is a relation between IR and UV 
renormalon residues which ensures the divergent term 
vanishes 



zl+sBiin + 3) = ~zlAi{n). 
This ensures that 

CxD 

J24[Mn) + B,{n)]^0. 

Another similar relation is [l3| 

Ao(n) = -Bo(n + 2). 



(93) 



(94) 



(95) 



We shall show that these relations are underwritten by 
continuity of the characteristic function in the skeleton 
expansion. At = A^ one finds the finite values, 

oc 

I?(^^(g2 = A2) = ^z„[^o(«)-So(«)] 

oo 



^z2[Ai(n) + Si(n)]lni 



0.679938 



(96) 



and 



/Cy^^(Q2^A2) = -^ln2 
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4T^(0'=A2) = --ln2. 



(97) 



5.1. QCD skeleton expansion 

In QED the insertion of chains of bubbles into a 
basic skeleton diagram produces a well-defined skele- 
ton expansion. In QCD by modifying the coefficient 
of In(-fcVQ^) in no(fc2) to involve the QCD beta- 
function coefficient h one can recast Eq.(71) in the 
form 




oo ,2 



fc2 \ fc2 



(98) 



which, defining t ~ k^/Q^, can be written as 



dt LUT,(t)a{tQ'^) 



(99) 



Here uj-D{t) is the characteristic function. It satisfies 
the normalization condition 



/"OC 

/ dtujv{t)^l 
Jo 



(100) 



which ensures that the leading coefficient of unity in 
the perturbative expansion of Eq.(31) is reproduced. 
ujT>it) and its first three derivatives are piecewise con- 
tinuous at ^ = 1 and the function divides into an IR 
and a UV part 

nl POO 

V^p^= / dtuj^^{t)a{tQ^)+ / dtuj^^a{tQ^mi) 
Jo J I 

We shall see that the IR k'^ < Q^) and UV (fc^ > Q2) 
components respectively reproduce the IR renormalon 
and UV renormalon contributions in the Borel plane. 
For > one encounters the Landau pole in the 
coupling in the first IR integral at t = A^/Q^, and the 
integral requires regulation (e.g. PV). This mirrors 
the IR renormalon ambiguities in the Borel integral 
of Eq.(74). For < one encounters the Landau 
pole in the coupling in the second UV integral, and 
this requires regulation. This mirrors UV renormalon 
ambiguities in the modified Borel integral of Eq.(90). 
uj-pit) can be derived by using classic QED work of 
Baker and Johnson on vacuum polarization. The vac- 
uum polarization function of Eq.(22) can be written 
as 



n(g') 



f 

Jo 



dt ujn{t)a{tQ^ 



(102) 



where the characterisitic function LUYi{t) is given by 
tE{t) t<l <-> IR 



ujn{t) 



(103) 



t> 1 ^ UV 



^(i) is given by 



.(t) . ^,[l-l-t+{---lnt)t 

+ ^i^^[L2(-0+ln<ln(l + <)]]. (104) 



where L2{x) is the dilogarithmic function 

ln(l - v) 



r 1 

L2{x) = - / dy- 
Jo 



y 



S(i) corresponds to the Bethe-Salpeter kernel for the 
scattering of light-by-light involving the diagrams in 
Fig. 6. Notice that by attaching the ends of the 
fermion bubble chain to the momentum k external 
propagators in Fig. 6 one reproduces the topology of 
the diagrams in Fig. 2. This computation of wn is 
a one-loop calculation, and as we shall see can be 
converted directly into the Borel plane renormalon 
structure for the D-function. This provides a much 
simpler route than the full two-loop calculation of 
Refs.mp. 

Changing from 11 to I? induces a transformation in 
ujii of 



n(Q') - Q'^n(Q2) = - 

UJjl{t) UJjl{t) + t^UJjl{t) 

dt 



One can write wn(i) as an expansion in powers of t 
and In t times powers of t 

. / oc oo \ 

\n=l n=2 ) 

and the conformal symmetry i ^ 7 between the UV 
and IR regions means that the UV part can be written 
in terms of the same coefficients 

. / 00 CXD \ 

= 5]e„i-"-ini^^;i-" (107) 



n=1 



The ^„ and ^„ are found to be 

4 (2-6ri2)(-l)^ 



Cn>l 



3 (n- l)2n2(n-|- 1)2' 

4 2(-l)" 

3 {n — l)n(n + 1) 



6 = 16 = 0. 



(108) 
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Figure 6: Light- by- light scattering diagrams, used to calculate ojn- 



For LOj) we have 

oo 

oo 
n=2 

oo 

oo 

+ lni^en(?^- (109) 

By making the change of variables z — —a{Q^){n + 
l)lni for i < 1 and z = a{Q^){n - l)lni for t > 1, 
one can transform the skeleton expansion result into 
the standard Borel representation of Eqs. (74,75). For 
q2 ^ a(Q^) < 0, and one obtains the modified 
Borel representation of Eq.(90) in which the upper 
limit of integration in z is — oo. For consistency one 
requires relations between the residues ^o,ii ^0,1, and 
the characteristic function coefRcients 



^„(l + n)+^„ 
n + 1 

^„(1 - n) - ^„ 



n- 1 

for the single pole residues and 
in{n + l) 



= -Bi{n + l)zn+i n>l 



Ai{n - l)zn-i n>2 (110) 



(n-h 1)2 

injn - 1) 
(n- 1)2 



= Boin + 1) + Bi{n + l)zn+i n>2 
= Aa{n - 1) - Ai{n- l)zn-i ?i(M2) 



for the double pole residues. These relations may be 
used to rewrite the series for LUvit) in terms of the 
residues ^o,ij -Bo.ii 

, oo 

oo 

- lni^(n + l)2[Bo(n + l) 

n=2 

+ Zn+iBi{n + l)]t" 



oo 

+ lnt^(n-l)'[Ao(n-l) 

n=2 

- Z„„l^l(71- l)]i-". 



(112) 



One can then show that continuity of to-Dit) and its 
first three derivatives at t = 1, and equivalently finite- 
ness of I?^(Q2) and its first three derivatives d/d\nQ 
at = A2 is underwritten by the following relations 
between the Ao,i and Bq^i residues 



J24[Mn) + B,in)]=0. 



(113) 



This is just Eq.(94) which guarantees finiteness of 
'D{Q'^) at (52 = 7^2. In addition we have more compli- 
cated relations which underwrite continuity and finite- 
ness of the derivatives. [H 

00 

Y.[2zl{A,{n) 

- Bi{n)) - zl{Ao{n) + Bo{n)] = . (114) 

OO 

J2[izt{A,{n) + B,{n)) 

- 2z3(^o(n)-Bo(n))] -0. (115) 

00 

Yy^zl{A^{n)-B^{n)) 

- 3ztiAo{n) + Boin))] = . (116) 



5.2. The NP component 

It is easy to show 1] that the ambiguous imaginary 
part in arising from IR renormalons for > 

A2 and UV renormalons for < can be written 
directly in terms of uj^ and Wp^, 



Imiv'j^^iQ^)] 



±- 



2Tr A2 



.IR 



2t: A2 



A2 



UV 



A2 



(Q^ > A2) 
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Continuity at = then follows from continu- 
ity of LUT>{t) at t = 1. In principle the real part of 
the OPE condensates are independent of the imagi- 
nary, but continuity and finiteness involve the set of 
relations between ^o.i and _Bo.i in Eqs. (113-116). Al- 
though not strictly necessary for continuity, continuity 
naturally follows if we write 



K± 



27ri 



dt ( UJT>{t) +t 



du)v{t) 
dt 



Here k is an overall real non-perturbative constant. If 
the PT component is PV regulated then one averages 
over the ± possibilities for contour routing, combining 
with V'j^^ one can then write down the overall result 
for 

Jo 



+ K UJT>{t)+t 



LOv{t)a{tQ^) 
dwv{t) 



dt 



e{A^ - tQ^j)] 19) 



The evolution is fixed by the non-perturbative 
constant k and by A. The evolution for the choices 
K = 0,1,-1 is shown in Fig. 7. We see that the 
QCD corrections to the parton model result correc- 
tions freeze smoothly to zero as ^ 0. 



6. Concluding remarks 

In these lectures we have shown how the large-order 
growth of perturbative coefficients in QED and QCD 
is dominated by renormalon contributions associated 
with diagrams such as those in Figs. 2, 4, in which 
chains of fermion bubbles are inserted inside a ba- 
sic skeleton diagram. We have focused on the vacuum 
-j^jjMlarization Adler function D(Q'^) as our main ex- 
arnple, with a cursory consideration of DIS sum rules 
as well. In QCD one is forced to use the so-called 
leading-^ approximation in which one inserts chains of 
effective bubbles in which the logarithm in the renor- 
malised no(A;^) has a coefficient of b, since correc- 
tions to the gluon propagator are gauge-dependent. 
In the approximation in which a single chain is in- 
serted one finds single and double pole IR and UV 
renormalon singularities in the Borel plane for Adler- 
D and for the DIS sum rules. Interestingly, these 
one-chain approximation results are finite and con- 
tinuous at = A^ where the coupling has a Landau 
pole, and this property is underwritten by relations 
between IR {k'^ < Q'^) and UV {k"^ > Q^) physics, 
where refers to the momentum flowing through the 
bubble chain. There are corresponding relations be- 
tween IR and UV renormalon residues. These con- 
nections are more transparent in the language of the 
one-chain skeleton expansion of Eq.(99) in which they 
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are seen to correspond to continuity of the character- 
istic function Lj{t) and its derivatives. Perturbative 
and non-perturbative ambiguities which contribute an 
imaginary part may also be written in terms of the 
characteristic function. Presumably the one-chain ap- 
proximation is far too crude to describe the real non- 
perturbative infra-red behaviour of QCD observables 
, but is interesting that it does have the properties of 
continuity and finiteness that must be possessed by 
the true all-orders result. It will be interesting in fu- 
ture studies to see how these properties arise if one 
includes higher numbers of chains, and also instanton 
effects. 
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